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Wc search for a light Higgs boson, AP , in the radiative decay T(35') ')Al\ AP t^t^ , 
— > e^v^r or — > \x^VpMT- The data sample contains 122 million T(3S') events recorded with 
the BABAR detector. We find no evidence for a narrow structure in the studied t'^t~ invariant 
mass region of 4.03 < < 10.10 GeV/c^ We exclude at the 90% confidence level (C.L.) a low 

mass Higgs decaying to r'^r~ with a product branching fraction B(T(3S') 7^°) x — > t"'"t~) 
> (1.5 — 16) X 10~' across the m^+^- range. We also set a 90% C.L. upper limit on the T"'"r~-decay 
of the r]h at B{r]h t+t~) < 8%. 
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PACS numbers: 13.20.Gd, 14.40.Gx, 14.80.Cp, 14.80.Mz, 12.60.Fr, 12.15.Ji 



In the standard model (SM) of particle physics [H, 
fundamental particles acquire mass via the Higgs mech- 
anism This mechanism requires the existence of at 
least one new particle called the Higgs boson. In the SM, 
there is only a single Higgs boson, with a mass of the or- 
der of the electroweak unification scale 100 GeV/c^). 
In the minimal supersymmetric standard model (MSSM) , 
additional Higgs doublets are required to give mass to the 
new particles (Sj. Moreover, in the next-to- minimal su- 
persymmetric standard model (NMSSM), an additional 
Higgs singlet field is introduced to solve the hierarchy 
problem A linear combination of this singlet with 
a Higgs doublet leads to a CP-odd Higgs state, A", 
whose mass need not be larger than 2mb, where mt is 
the &-quark mass [1, [1]. It is ideal to search for this 
state in T ^ decays @. The branching fraction 
B{r{'3S) jA") depends on the NMSSM parameters, 
but a value as large as 10~^ is plausible for reason- 
able parameters In the mass range where the decay 
A° — > r+T~ is kinematically accessible, this mode is ex- 
pected to dominate. Constraints on the invisible and 
dimuon decays of the A" have recently been obtained. 

The current best limit on the product of branching 
fractions B{T{IS) -fA°) x B{A° r+r") is given 
by the CLEO Collaboration foj based on a data sample 
of 21.5 million T{1S) candidates. The CLEO 90% C.L. 
limits cover the range 2mr < mj^a < 9.5 GeV/c^ (toi- is 
the r-lepton mass [10]) and vary between 1 x 1Q~^ and 
48 X 1Q~^. A recent DO search for a neutral pseudoscalar 
Higgs boson in a similar mass range showed no significant 
signal fll,]. 

In this Letter, we study the decays T{3S) 7r+T^, 
where the search for A^ is extended for a wider mass 
range w.r.t the T{1S) jt^t^. We scan for peaks 
in the distribution of the photon energy, E^^ corre- 
sponding to peaks in the rr invariant mass m'^+^- — 
7/135 ~ '^mssE.y, where 17133 is the T{3S) mass and is 
measured in the T{3S) rest frame (center-of-mass (CM) 
frame) . We quote branching fraction values in the region 
4.03 < TOt-+t-- < 10.10 GeV/c^, but we exclude from our 
search the region 9.52 < m^+T- < 9.61 GeV/c^, because 
of the irreducible background of photons produced in the 
decay chain r(35) ^ iXbji^P), Xbji^P) ^ lT{lS), 
where J = 0, 1, 2. In addition, we set an upper limit on 
B{r]b -> r+T"). 

The data were collected with the BaBAR detector 12 1 
at the PEP-II asymmetric-energy e~^e~ storage rings at 
the SLAC National Accelerator Laboratory, operating at 
the T{3S) resonance. We use a data sample of 122 million 
T(35) events, corresponding to an integrated luminosity 
of 28 fb~^. We also use data samples of 2.6 fb^^ recorded 
30 MeV below the T(3S') (0FF3S), 79 fb-^ at the T{AS) 
(0N4S), and 8 fb"! 40 MeV below the r(45) resonance 



(0FF4S) to study the background and to optimize the se- 
lection criteria. These data samples were taken with the 
same detector configurations. Monte Carlo (MC) event 
samples based on Geant4 [l^] simulation of the detec- 
tor are used to optimize selection criteria and evaluate 
efficiencies. 



We select events in which both r-leptons decay lep- 
tonically, r+ — > e'^i'eVr or r"*" —>■ ^'^Vj^t (denoted in 
the following as r ^ e, or t ^ fj.) [ij]. Events are 
required to contain at least one photon with E^> 100 
MeV, and exactly two charged tracks. We allow up to 
nine additional photons with energies below 100 MeV in 
the CM frame. Photons are reconstructed from local- 
ized deposits of energy in the electromagnetic calorime- 
ter, which have energies larger than 50 MeV in the labo- 
ratory frame and which are not associated with a charged 
track. Both charged tracks are required to be identified 
as leptons (e or /i) . After this selection the residual back- 



ground is mostly due to 



jT~^T and higher order 



QED processes, including two-photon reactions such as 
e^e^ — > e+e^e+e^ and e~^e~ — > e^"e^/i+/i^ with smaller 
contributions from other T{3S) decays and e^e^ — > qq 
{q = u,d,s,c). 

To reduce this residual background, we exploit a set 
of eight discriminating variables: the total CM energy 
(^-totai) calculated from the two leptons and the most 
energetic photon; the squared missing mass (rnf^^^s) 
tained from the missing four-momentum, which is the dif- 
ference between the final and initial state momenta; the 
aplanarity (Api), which is the cosine of the angle between 
the photon and the plane of the leptons; the largest cosine 
between the photon and one of the tracks (cos 0-y_track); 
the cosine of the polar angle of the highest-momentum 
track (cos^track); the transverse momentum of the event 
(pt) calculated in the CM frame; the cosine of the polar 
angle of the missing momentum vector (cos^miss); and 
the cosine of the opening angle between the tracks in the 
photon recoil frame (cos^opon)- The final selection cri- 
teria on these variables are obtained by maximizing the 
quantity S/\/B, where S (B) stands for the expected 
number of signal (background) events. Numbers of sig- 
nal events are obtained from MC samples, while back- 
ground yields are obtained from the 0FF3S, 0N4S, and 
0FF4S datasets. Since the background varies as a func- 
tion of the photon energy, we optimize the selection cri- 
teria in five E^ regions: (Si) 0.2 < E^ < 0.5 GeV, (6*2) 
0.5 < E^ < 2.0 GeV, (5*3) 1.5 < E^ < 2.5 GeV, (5*4) 
2.5 < E^ < 3.5 GeV, and {S5) 3.0 < E^ < 5.0 GeV. 
The overlaps between the Ej regions reduce the discon- 
tinuity in the efficiency at the boundaries. The dominant 
irreducible background is due to e~^e^ — > 7t+t^. The 
highest level of background contaminations is observed at 
low E^ values. Among the different final states, the back- 
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ground is largest in rr ee and smallest in rr e/i. 

The photon energy resolution degrades as a function of 
E^, from 8 MeV at E"^- 0.2 GeV to 55 MeV at E^'^ 4.5 
GeV. The selection efficiency is calculated using MC 
events. The efficiency in the rr ee, rr — > e/z, and 
rr — > /i/i modes varies as a function of E^ between 10- 
14%, 22-26%, and 12-20%, respectively. The MC sam- 
ples are generated with angular decay distributions ex- 
pected for a CP-odd Higgs boson; similar efficiencies are 
obtained for CP-even states. 

We search for an excess in a narrow region in the E^ 
spectrum since any peak in the recoil mass {rrirr), indi- 
cating the presence of a new particle decaying in r-pairs, 
translates to a peak in the E^ distribution. We describe 
the E^ distribution as a smooth background spectrum 
and a narrow enhancement of known width, but unknown 
position and event yield. We perform a binned maximum 
likelihood fit simultaneously to the rr — > ee, e/i, and /i/i 
samples. 

The fit is performed in two steps. First, we assume 
there is no signal and fit the background function. The- 
oretical motivations [15] inspired the choice of the back- 
ground function shape, / — [p{l — xY / Ei^ + s / E^) ■ 

/3(x)-(3-/32(a;)), where /3(a;) = ^ I - Ami/ {mlg{l - x)) , 
X = 2E^/'m3S- For each rr-decay mode, a different set 
of the parameters p, q, r, s is used. These parameters are 
allowed to vary. 

The events r(35) ^ 7X6j(2P), Xbj{2P) ^ ir{nS), 
and T{nS) t^t^ (J = 0,1,2; n — 1,2) are expected 
to peak in E^ when the photon from X6j(2P) ^T{nS) 
is misidentified as the radiative photon from the T(iS) 
decay. Each of the peaks in the photon spectrum due 
to the \tj(2P) ^T{1S) transitions is described by a 
Crystal Ball [la| (CB) function. The mean values for 
the Yf,n(2P) and Xf,i(2P) CB functions are fixed to the 
PDG \ld\ and the width values are fixed to the MC 
resolution, while the mean and width for Xf,2(2P) are 
free. The power law and the transition point for all 
CB functions used in the analysis are fixed to the val- 
ues obtained in MC. The event yields for the Xbj{2P) 
background for each of the three rr data samples are 
related via their relative efficiencies, which are func- 
tions of E^. To account for the contributions from 
XbjC^P) — > 7^^(25'), a fourth CB function is added, for 
which the mean, width, and the relative normalization 
are free. The fitted mean and width obtained for this 
peak are 234 ± 2 MeV and 13.3 ± 2.7 MeV (statistical 
uncertainties only), respectively. The number of events 
from the Xf,j(2P) ^T{nS) [n — 1,2) contamination 
are common between the different rr-decay modes, and 
divided between these modes according to the efficiency 
sum, = -I- 2ee^ -I- e^^, where egg, Cg^, and e^^ are 
the efficiencies as a function of E^ in the decay modes 
rr — > ee, e/i, and /i/x, respectively. An example of the 
fits to the E^ distributions in the different r+r~-decay 



modes, obtained with the selection criteria and fitted 
in the region 0.2 < Ej < 2.0 GeV, are shown in Fig. [1] 
Satisfactory fits are obtained. 
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FIG. 1: (a), (c), (e): Ej distributions for the differ- 
ent rr-decay modes. Filled circles show the data; dot- 
ted lines represent contributions from T{3S) — » 7X6j(2P), 
Xbj(2P) — > ^T{2S); dotted-dashed lines show contributions 
from r{SS) 7X(,j(2P), Xi>j(2P) ^ 7^(15'); and solid lines 
show the total background function. For each rr-decay mode, 
the difference between the background function and the data 
divided by the uncertainty in the data is shown ((b), (d), (f)). 

In the second step of the fit procedure, we search for 
the signal T{3S) 7^°, A" t+t~. We assume the 
A° has negligible width [TTj , and parameterize the signal 
distribution with a CB function. The search for such a 
signal is performed by scanning for peaks in the E-y distri- 
butions in steps that are equal to half the photon-energy 
resolution at any chosen value of Ej. In total, 307 scan 
points are examined. The mean of the signal function 
is fixed to the photon energy at the i"^ scan point (El^). 
The signal width is fixed to the value of the photon en- 
ergy resolution obtained from the MC simulation. The 
contribution from each rr-decay mode to the total num- 
ber of Higgs candidates is proportional to the fractional 
efficiency for a particular mode. The background shape 
parameters (including the XbJ parameters) are all fixed 
to the values determined in the first step of the fit, with 
the exception of p and s, to allow free background nor- 
malization. The number of free parameters in each fit is 
seven (pee, Pep,, Ppp, See, Sep, Spp, and iVsig), where the 
subscripts indicate the final state of the rr-decay modes. 
When the scan is performed in the regions 5*3, 5*4, and 
S5, the parameters See, Sep, and s^^ are fixed to zero. 
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For each scan point, the yield, 
cal uncertainty, a{Nsig) are obtained from the fit 



A^sig, and its statisti- 



The 



yield significance from the data, Nsig/a{Nsig) is shown 
in Fig. [21 and overlaid with a standard normal distri- 
bution. The data points are consistent with the normal 
distribution, and therefore no significant evidence for any 
unknown narrow structure is observed in the scan. 
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FIG. 2: Nsig/a'{Nsig) as obtained from the scanning proce- 
dure. Only statistical uncertainties are included. The curve 
shows the standard normal distribution with a normalization 
factor of 307. 



Product branching fractions are determined from the 
signal yields at each scan point, correcting for a fit bias 
described below. The results are shown in Fig. [3fa). 
These results show no evidence for a narrow resonance 
in the mass range under study. Bayesian upper limits 
on the product of branching fractions, computed with a 
uniform prior at 90% C.L., are shown in Fig. [3]Jb). The 
solid line shows the limits obtained with the total uncer- 
tainties (statistical and systematic added in quadrature) 
while the dashed line shows the limits with statistical 
uncertainties only. 

We measure the branching fraction Bijjb — > r+r^) — 
(-0.1 ±4.2 ±2.3) % at m^+^= 9.389 GeV/c^, using the 
B{T{3S) -fr]b) from Ref. Therefore, the 90% C.L. 
upper hmit on B{i]i, t^t~) is 8 (7)%, considering all 
(statistical only) errors and accounting for the expected 
10 MeV width of the rib- We note that the limit and 
branching fraction are insensitive to the rj^ width within 
the expected 5-20 MeV range [ll|. 

We account for systematic uncertainties due to track- 
ing (2%), lepton identification (1.2-2.6%, depending on 
the TT-decay mode), photon reconstruction efficiency 
(4%), and the number of T{2,S) (1%). In the scan pro- 
cedure, the parameters of the background shape and of 
the Xbj{2P) states are fixed. To estimate the systematic 
uncertainty related to these parameters, each parameter 
is varied by its estimated statistical uncertainty deter- 
mined in the first step of the fit. The scan procedure is 
repeated for each parameter change. When calculating 
the systematic uncertainties from this source, the cor- 
relations between the various parameters are taken into 




FIG. 3: (a) Product branching fractions as a function of the 
Higgs mass. For each point, both the statistical uncertainty 
(from the central value to the horizontal bar) and the total 
uncertainty (statistical and systematic added in quadrature) 
are shown (from the central value to the end of the error 
bar). In (b), the corresponding 90% C.L. upper limits on 
the product of the branching fractions versus the Higgs mass 
values are shown, with total uncertainty (solid line) and sta- 
tistical uncertainty only (dashed line). The shaded vertical 
region represents the excluded mass range corresponding to 
the Xf)j(2-P) -yTilS) states. 



account. The ratio between the total systematic uncer- 
tainties due to the background shape and the statistical 
uncertainties varies between 12% and 170%. The largest 
systematic variations occur for larger values of m^+^~, 
and are due to the uncertainty in the q^fj, parameter for 
TT e/i. The fit bias and its uncertainty are determined 
by applying the fit procedure to a large number of MC 
experiments. Each MC sample contains a known number 
of signal events, while background events are generated 
according to the background shape. The event yield, re- 
turned by the fit, is a linear function of the number of in- 
put events. The event yield in the data is corrected using 
this function. The difference between the corrected and 
uncorrected event yield is (conservatively) considered as 
the systematic uncertainty due to the fit bias, which is 
typically small (few percent) but can be as large as 30% 
of the statistical uncertainty at high to^+^- . The system- 
atic uncertainty associated with the choice of the signal 
shape function is determined by varying the values of the 
parameters in the signal CB function; the width and the 
power law are varied (multiplicativcly) by 30% and 38%, 
respectively; the transition point is varied (additively) by 
36%. The associated systematic contribution is typically 
small (few percent) but is as large as 50% of the sta- 
tistical uncertainty at large m^+^- . Finally, we include 
a systematic uncertainty of 0.6% to account for the sys - 
tematic uncertainty due to the r branching fractions [1Q| . 
The dominant systematic uncertainties are due to the 
background-shape parameters, which are obtained from 
fitting the same data sample. Thus, we conclude that the 
main systematic uncertainties are primarily statistical in 
nature. 

In summary, we have performed a search for a light 
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Higgs boson in the radiative decays T{3S) jt^t~, 
where — > e^VgVr or r"*" fX^v^Vr, using a data sam- 
ple of 122 miUion T(3S') events. Our search covers the 
mass range 4.03 < m^+T-- < 10.10 GeV/c^, excluding 
9.52 < m^+^- < 9.61 GeV/c^ to veto the Xbji'^P) with 
Xf,j(2P) jTllS). No evidence for a signature of light 
Higgs boson decays to r pairs is observed. In this mass 
interval, the upper limits on the product branching frac- 
tion B{r{3S) 7^°) X B{A° T+T-) vary between 
(1.5 - 16) X 10-5 at 90% C.L. 
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